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SOURCE DISTRIBUTION METHOD FOR UNSTEADY ONE -DIMENSIONAL FLOWS WITH 
SMALL MASS, MOMENTUM, AND HEAT ADDITION AND SMALL AREA VARIATION 

By Harold Mirels 


SUMMARY 

A source distribution method is presented for obtaining flow per- 
turbations due to small unsteady area variations, mass, momentum, and 
heat additions in a basic uniform (or piecewise uniform) one -dimensional 
flow. First, the perturbations due to an elemental area variation, mass, 
momentum, and heat addition are found. The general solution is then 
represented by a spatial and temporal distribution of these elemental 
(source) solutions. Emphasis is placed on discussing the physical nature 
of the flow phenomena. 

The method is illustrated by several examples. These include the 
determination of perturbations in basic flows consisting of (l) a shock 
propagating through a nonuniform tube, (2) a constant- velocity piston 
driving a shock, (3) ideal shock-tube flows, and (4) deflagrations in- 
itiated at a closed end. The method is particularly applicable for find- 
ing the perturbations due to relatively thin wall boundary layers . 


INTRODUCTION 

Numerous studies have been made of unsteady one -dimensional flows 
through ducts with area variation, body forces, mass, and heat addition. 
The most general method of handling such problems is by ^the method of 
characteristics which, for these problems, employs the "Riemann variables. 
A full account of this method is given in reference 1. 

When the area variation, body forces, and so forth introduce only 
small disturbances into an otherwise uniform (or piecewise uniform) basic 
flow, the equations of motion can be readily linearized. As is usual 
with linear problems, many different methods and viewpoints can then be 
used to obtain solutions. One approach is to find first the perturba- 
tions associated with an elemental area change, body force, mass, and 
heat addition. The general linearized solution of an unsteady one- 
dimensional flow can then be represented by a spatial and temporal 



2 


distribution, of such elemental solutions . rhis approach may be termed 
a "source distribution method." The advantage of this viewpoint over 
more formal methods of solution is that it aften gives a better physical 
insight into the flow processes. A discussion of the source distribu- 
tion concept , as applied to one -dimensional unsteady flows, does not 
appear to exist explicitly in the current literature. Hence, such a 
discussion is presented herein. To emphasise the physical nature of 
these flows, the solutions are first deduced from a relatively simple 
flow problem. The solutions are then again deduced by formal' manipula- 
tion of the equations of motion as they appear in the method of charac- 
teristics. The use of these solutions is then illustrated by several 
examples . 


The source distribution method has been previously applied, by the 
author, to find nonuniformities in a shock uube due to the unsteady 
boundary layer along the shock tube wall (refs. 2 and 3). The present 
report may be considered as an extension and elaboration of the method 
presented in those reports . 


ANALYSIS 

A uniform basic flow is assumed to be £ lightly disturbed by small 
area changes, body forces, mass, and heat additions, which induce un- 
steady one -dimensional perturbations. The j erturbations due to an ele- 
mental area change, body force, and so forth are first deduced from a 
simple flow problem. The general solution is then expressed as the 
superposition of these elemental flows. Finally, the general solution 
is again obtained by a formal linearization of the equations of motion 
as they appear in the method of characteristics. 


Equations of Motion 

Consider a uniform basic flow through a tube of constant area A. 
The uniform fluid properties are denoted by the symbols p, p, T, 
u -> ' ‘ • • (Symbols are defined in appendix \.) Let perturbations from 
those uniform values be denoted by the prefix A so that the net pres- 
sure at a point is p + Ap, the net density is p + Ap , .... The per- 
turbations are generally functions of (x,t) 3 o that Ap = Ap(x,t), 

Ap = Ap(x,t) , and so forth. By assuming that the perturbations are due 
to small area variation, body forces, mass, and heat additions, the 
equations of motion are: 

Continuity: 


^Yp 

+ 


dAp cAu 

"T -f p -T* 

dx Ox 


:= P 


(la) 



Momentum : 
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Energy: 


State : 


cAu dAu 1 cAp f 

■5r + u ^r + p^T'p 


dAs dAs q 

"5t *" u "ST = "pT 


Ap Ap AT 
p " p ' T 


(lb) 

(lc) 

(l d) 


where 

, . P /^AA. dAA.\ 

p(x,t) = m - ^ 


m(x,t) = mass addition, per unit volume, per unit time 


AA.(x,t) n perturbation of cross-sectional area 

f(x,t) = body force, per unit volume, acting in +x direction 


q(x,t) 2 heat addition, per unit volume, per unit time 


As(x,t) z c 




entropy perturbation 


The quantities p, f, and q may be referred to as volumetric sources. 

They are, respectively, sources of mass^", momentum, and heat. The term 
"volumetric" stems from the fact that they are defined on a "per unit 
volume” basis. 


For a given area variation AA.(x,t) , the average velocity (normal 
to the wall) of a fluid particle at the wall at any point (x,t) is given 
dy 

1 ( cAA , 2 aa\ 

v “ - t [sr + u 


•^The quantity p may be considered as the equivalent volumetric 
mass source in a tube of constant area, which induces the same pertur- 
bations as a prescribed mass addition m and area variation AA. 
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where l is the perimeter corresponding to the cross-sectional area A. 
The normal velocity is positive when directed inward. Then, for the case 
of area variation (but no external mass addition). 


P = 


Pi 

A 


v 


( 2 ) 


This expression is convenient for finding perturbations due to thin un- 
steady boundary layers along a constant-area tube. Here, v is the ver- 
tical velocity at the outer edge of the boundary layer, as computed from 
boundary-layer theory. This approach is used in references 2 and 3 to 
nonuniformities in shock tubes due to the wall boundary layer. 


Elemental Source Solutions 

The quantities p, f, and q were referred to as volumetric sources 
since they were defined on a per unit volume basis. Equivalent quanti- 
ties, defined on a per unit cross-sectional area basis, may be referred 
to as "planar sources." For example, plana.* mass, momentum, and heat 
sources at x = 0 can be defined by the relations 


p = lim 

/ H dx 

f = lim 

f fdx 

q = lim 

€->0 ’ 

- € 

€-+0 * 

-€ 

e-MD 

p-*oo 


f+oo 


q+oo 


/: 


q dx 


(3) 


where 


p = equivalent mass addition at x = 0, per unit A, per unit time 
f = body force at x = 0, per unit A, in direction 
q — heat addition at x = 0, per unit A, per unit time 

In general, p, f, and q can be functions cf time. In the present sec- 
tion, a simple flow containing planar sources will be considered. The 
solution of this flow will then be used to ceduce the flow field induced 
by an elemental source. 


Assume that planar mass, momentum, and heat sources are placed at 
x = 0 at time t = 0 (in an otherwise uniform flow) and that they re- 
main at constant strength thereafter. The resulting perturbed flow field 
is indicated in figure 1 (for the case where the basic uniform flow is 
subsonic). In tigure l(b), line a Is a downstream propagating acoustic 
wave moving with velocity u + a, line b is a contact surface between 
two regions of different entropy, moving witn velocity u (neglecting 
perturbation velocities), and line c is an uostream propagating acoustic 



5 


wave moving with velocity u - a. The lines a, b, c , and x = 0 sub- 
divide the flow into four regions , 1, 2, 3, and 4. Region 1 is the 
original undisturbed uniform flow. The perturbations in regions 2, 3, 
and 4 are independent of x and t. 

The magnitude of the perturbations in the various regions can be 
found by considering the jump in the perturbation quantities across 
lines a, b, c, and x = 0. The case M< 1 is treated here, while the 
case M > 1 is treated in appendix B. 

Line a. - Line a is a downstream propagating acoustic wave. The 
acoustic relations then give 


Ap 2 = paAu 2 (4a) 

Ap 2 = Au 2 (4b) 


Line b. - The fluid velocity and pressure are continuous across a 
contact surface so that 


Au 3 = Au 2 (5a) 

Apg = Ap 2 (5b) 

Line c. - Line c is an upstream propagating acoustic wave. The 
acoustic relations give 

Ap 4 = -paAu^ (6a) 

Ap 4 = - •“ Au 4 (6b) 

Line x = 0. - Integrating equations (la) to (lc) across x = 0 
and noting equations (3) give, respectively, 

u(Ap 3 - Ap 4 ) + p(Au 3 - Au 4 ) = P (7a) 

-i £* 

u(Au 3 - Au 4 ) + ± (Ap 3 - Ap 4 ) = — (7b) 


u 




q. 

c v pT 


Equation (7c) follows from As 4 = 0, As 3 = c v [(Ap 3 /p 3 ) - r(Ap 3 /p 3 )]. 


(7c) 
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Equations (4) to (7) are nine equations in nine unknowns. The solu- 
tion can be expressed as 


Ap 2 _ 

r 

(- . q , t \ 

(8a) 

P 

2(1 + M)pa 

V C P T a ) 

Ap 4 

r 

|ch 

1 & 

(8b) 

p 

2(1 - M)pa 


As 3 

c v 

q 

c v puT 


(8c) 


AH other perturbations are then found directly from equations (4) to (7) 

Now , assume that equal and opposite ( v ith respect to the previous 
example) planar mass, momentum, and heat sources are placed at x = 0 
at time t = At (in the same otherwise uniform flow) and that they re- 
main at constant strength thereafter. The corresponding t-x diagram 
is shown in sketch (a). The lines representing the acoustic waves and 



the contact surface are displaced upward by an amount At from the cor- 
responding lines in figure l(b). The perturbations in regions 2', 3', 
and 4 are equal and opposite to those of the corresponding regions in 
figure 1(b). 

If the flow fields in figure l(b) and sketch (a) are linearly super' 
posed, the resulting flow (sketch (b)) has zero perturbations everywhere 
except between lines a, a' j b, b 1 j and c, c ' (since the perturbations in 
regions 2, 3, and 4 are equal and opposite to those in regions 2', 3', 
and 4'). By using a superscript ( ) + to denote perturbations due to a 
downstream propagating acoustic wave, the pressure perturbations between 
lines a and a' can be written (from eq. (8a)) 
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/ 



(b) 



p 


r 

2(1 + M)pa 



(9a) 


Also, for this region (from eqs. (4)), (Au + /a) = (Ap + /p) = (Ap + /Tp) ■ By 
denoting the perturbations due to an upstream propagating acoustic wave 
by the superscript ( )~ , the pressure perturbation between lines c and 
c 1 can be written (from eq. (8b)) 


ap~ r (- , _ji_ t\ 

p 2(1 - M)pa V c p T a ) 


(9b) 


Also, for this region (from eqs. (6)), (-Au"/a) = (Ap - /p) = (Ap"/rp) • 
Between lines b and b', the velocity and pressure perturbations are zero, 
while the entropy perturbation is given by (from eq. (8c)) 


As _ q 
c v ~ c v puT 


(9c) 


The corresponding density perturbation is (Ap/p) = -(As/c p ). 

A similar development is presented in appendix B for the M > 1 
case. The results for Ap + , Ap" , and As are the same as for the M < 1 
case (e.g., eqs. (9)) except that the expression for Ap” contains the 
term (M - l) for the M > 1 case. Use of |l - M| in the demonimator 
of equation (9b) makes equations (8) applicable for both M < 1 and 
M > 1. 


Now, consider sketch (b) in the limit as At -* 0- It is then clear 
that a planar source at x = 0, t = 0 instantaneously generates a 
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downstream propagating acoustic wave, an upstream propagating acoustic 
wave, and an entropy wave which propagates, along the lines x = (u + a)t, 
x = (u - a)t, and x = ut, respectively. The strength of these waves is 
given by equations (9). The perturbation:: are zero everywhere except on 
these lines . ~" 

The results of the previous paragraph may be readily generalized to 
the case where the planar source is placet, at the arbitrary point x = 
t = r (see sketches (c) and (d)): 


M < I M > I 



The resulting perturbations and the lines along which these perturbations 
propagate are then 




(10a) 


(10b) 


along x=!+(u-a)(t-r), and 
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As _ q 

c y “ c y puT | 

= ~~ (Ap = Au = 0) J (lOc ) 

along x = | + u(t - t). 

The perturbations are zero everywhere except on these lines. The lines 
along which the disturbances propagate are termed "characteristic 1 ’ 
lines. The line along which the entropy wave propagates may also be 
referred to as a "particle path" line since the entropy wave is convected 
by the stream. 

Equations (10) are the elemental source solutions_which we had set 
out to obtain. From these equations it is seen that \i and q generate 
upstream and downstream propagating acoustic waves which have the same 
sign, whereas f generates upstream and downstream propagating acoustic 
waves which have opposite signs. The ratio of the downstream and up- 
stream acoustic waves is Ap + /Ap” = |l - M|/(l + M) for mass or heat 
addition and is Ap + /AP~ = - |l - M|/(l + M) for a body force. Thus, the 
upstream propagating acoustic wave is stronger than the downstream wave 
(considering p, q, and f separately) except for M = 0 and M 
As M approaches 1, the value of Ap" becomes very large (violating 
the assumption of small perturbations) and equations (10) become invalid. 


General Solution 

The perturbations due to an arbitrary spatial and temporal distribu- 
tion of volumetric mass, momentum, and heat sources can be found by the 
linear superposition of the elemental solutions of equations (10) . Again, 
use the coincident coordinate systems (£,r) and (x,t) where (£,t) defines 
the source location and (x,t) is the point at which the perturbations are 
to be found. Note that the volumetric mass source between £ - (d£/2) 
and £ + (d|/2) is equivalent to a planar mass source of strength 
p(£,t) = p(|,T)d£, and so forth. Then, by linear superposition, the net 
perturbation at (x,t) due to an arbitrary volumetric source distribution 
is given by 


Ap(x,t) _ Ap + (x,t) + Ap~(x,t) 

P P P 

. Au(x,t) _ Ap + (x,t) Ap~(x,t) 


(11a) 


a 


P 


P 


(in) 
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Ao(x,t) _ Ap(x,t) £ s(x,t) 
P TP ‘ c p 

AT(x,t) _ Ap(x,t) Ap(x,t) 
T P “ P 


(11c) 

(lid) 


where 


Ap + (x,t) 


T(T+l4jpa 




T=t - 


xzj 

u+a| 


id£ (lie) 


Ap~(x,t) 

P 




2(1 - M)pa 


As (x,t) 




T=t - 


X-| 

u-a I 


ks (nf) 


c v puT 


d| 


(Hg) 


The integrations in equations (lie) to (llg) add the contributions of 
all the sources which contribute to Ap + , A)", and As at (x,t). These 

v M t 

sources lie along the characteristic lines t = t an( j 

x _ | u ± a 

t = t - as indicated in sketches (e) and (f ) : 


M < I 


M> I 



The upper limit in equation (llf) is (-n») or (-<») depending on whether 
M < 1 or M > 1. (By iising this convention, d£/(l - M) is always posi- 
tive, and the absolute value sign is not needed for 1 - M. ) 
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Substitution of equations (ll) Into equations (l) verifies that they 
are indeed, the general solution of equations (l) • Equations (ll) could 
have been deduced, formally, from equations (l) • (This is, in fact, done 
in the next section.) The present development was undertaken so as to 
hring out more clearly the physical nature of the solution. Equations 
(lie) to (llg) assume that all perturbations are due only to the speci- 
fied p, q, and f distribution, so that there are no extraneous waves 
propagating in the tube. Otherwise, arbitrary functions of x - (u + a)t, 
x - (u - a)t, and x - ut (i.e., homogeneous solutions of eqs. (l)) would 
have to be added to the right sides of equations (lie) to (llg), respec- 
tively. The quantities p, q, and f must approach zero sufficiently 
fast, as | approaches l<», to make the integrals converge. As a result, 
the perturbations at x = j?° must be zero for all problems where equa- 
tions (11) apply. 


If the integrations in equations (lie) to (llg) are made with respect 
to i instead of £, the integrals take on the form (since 
dx/dl = l/(ui a) for t = t - (x - t)/{u± a) and dx/d£ = l/u for 
t = t - (x - |)/u) 



(12a) 


(12b) 



(12c) 


Equations (12) are somewhat simpler than equations (lie) to (llg). The 
integrations always proceed in the +x direction. Note that the coef- 
ficient 1 - M does not appear in equation (12b) . Thus, Ap does not 
become infinite as M - 1 (provided p, q, and f approach zero suffi- 
ciently fast, as T -> —co ) y and equations (12) are generally applicable 
at M = 1 as well as M ^ l. 2 (Eqs. (ll) are also applicable at M = 1 
provided d^/(l — M) is treated as an indeterminate form. ) 

%his contrasts with the earlier result that Ap °° as M ->■ 1 
for a planar source (eq. (lOb)). This difference between the perturba- 
tions due to a planar source and a volumetric source is typical of source 
distributions in fluid flow problems. The intensity of the singularity 
induced by a source in a fluid flow field decreases when a point source 
is replaced by a surface distribution of sources and when a surface dis- 
tribution is replaced by a volume distribution. 
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Alternate Derivation of General Solution 


The general solution obtained in the previous section is rederived 
herein by a formal manipulation of the equations of motion as they appear 
in the method of characteristics (ref. l). The latter employs the 
Riemann variables P and Q as dependent variables. The Riemann vari- 
ables are then related to the dependent variables Ap + and Ap“ . 


To solve one -dimensional flows of a perfect gas (with area change, 
body forces, mass, and heat addition) by the method of characteristics, 
the equations of motion can be written in the form (from ea. (ill d 9) 
of ref. 1) ' ' 


where 


8 + P a 

„ + -1- + 1 + P 6 + S l 

5t 2p 

C P T a C p^ r " 5t J 

a 

M + 1 - f + P 6 - S l 

5t 2p 

c p T a c pfr - 1) st J 


Ds 

Dt “ pT 




a( ) 

dt 


+ (u ± a i 



(13a) 

(13b) 

(13c) 

(13d) 

(I3e) 

(l3f ) 


P( ) 3 ( ) , 

Dt “ “ot - 



(I3g) 


The quantities P and Q are the Riemann variables. (These equal one- 
half the values as defined in ref. 1. ) Equations (13) define the varia- 
tion of P, Q, and s in the characteristic directions dx/dt = u + a, 
u - a, and u, respectively. A numerical integration of these equations 
(together with the equations of state) can then be obtained by proceed- 
ing along the characteristic directions, in small increments, as dis- 
cussed in reference 1. 


H Pj 1 > and. q are small, the equations can be linearized. Again, 
by letting the prefix A represent the departure of a flow variable from' 
the basic uniform flow, equations (13) become 
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A 

AP 

As 

r 

St 

a 

i — i 

i 

ft 

o 

oo 

“ 2p 

s_ 

AQ _ 

As 

r 

St 

a 

i 

f — % 

i — t 

i 

ft 

o 

(XI 

2p 



BAs _ q 
Dt “ pT 


(14a) 

(14b) 

(14c) 


Integrating equations (14a) and (14b) along the characteristic directions 
gives (for AQ - AP = As = 0 at t = -oo) 


r 


As 

2c p (r - i) 


]-h 


. aq 

a 


As 

a 2c p (T - 1) 


X 

2p 


J h(5,t) ^ 


^dT 

|=x-(u+a) (t-r) 

(l5a) 
dT 


=x- (u-a) (t-T ) | 

(15b) 


The integral of equation (l4c) is the sane as equations (llg) or (12c). 
Comparison of equations (15) with equations (12a) and (12b) shows 


r 

r 


ap 

a 

AQ 

a 


As 

2c p (r - i)_ 

n 

As 

2c (r - i) 

p _ 


II 

+ 

(16a) 

11 

1 

(16b) 


Thus, equations (12) (and eqs . (lie) to (llg)) could have been obtained 
by a formal integration of equations (13a) to (l3c). Also, the elemental 
source solutions can be deduced from equations (15). However, it is felt 
that the previous derivation of these equations brings out the physical 
nature of the flow. 


Equations (16) relate the perturbation of the Riemann variables 
AP and AQ to the variables of the previous section (namely Ap + , Ap - , 
and As). This relation can be deduced directly by a perturbation of 
the Riemann variables. Thus, perturbing equations (13d) and (l3e) gives 
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r Jap 

a | AQ 



Ap + rAs 

P 2 Cp (r - 1) 

Ap~ rAs 

p 2c P (r - i) 



+ 


+ 


rAs 

2c p (r - i) 
rAs 

2c p (r - l) 


(17) 


(18) 


Equations (17) agree with equations (16), as expected. Equations (17) 
and (18) provide a physical interpretation for AP and AQ within the 
limitation of linearized theory. These equations show that for isen- 
tropic flow (i.e. , As = 0) the local values of rAP/a and rAQ/a ex- 
actly equal Ap + /p and Ap"/p, respectively. Or, alternately, the local 
values of AP and AQ equal Au + and -Au", respectively. For As ^ 0, 
the relation between AP, AQ and Ap + , Ap" (or Au + , Au") also depends 
on the local value of As and, therefore, depends on the thermal history 
of the fluid element at the section under consideration. 


A linearized method of characteristics, employing ZIP, AQ, and As 
as the dependent variables, was developed in reference 4 for finding non- 
uniformities in a shock tube (due to the wall boundary layer). The method 
of reference 4 can be compared with the prasent method as follows. Ref- 
erence 4, in effect, integrates equations (13a) to (13c) in the appro- 
priate characteristic directions. Since s appears on the right side of 
equations (l3a) and (l3b) , the solution for AP and AQ is coupled with 
the energy equation (eq. (l3c)) for nonise: itropic flows, and all three 
equations must be solved simultaneously. In the present method, Ap + , 

Ap", and As are the dependent variables. However, the solution for 
Ap + and Ap" is not coupled with the energy equation and can be found 
without a knowledge of As. Thus, the solution for Ap + and Ap" is 
somewhat simpler than the corresponding solution for AP and AQ. In 
addition, Ap + and Ap" have a simple physical interpretation, which is 
not the case for AP and AQ. For isentropic flow, the two methods be- 
come identical (except for notation) . 
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APPLICATIONS 

The applications of equations (ll) and (12) are illustrated by con- 
sidering several one-dimensional unsteady flow problems. In particular, 
three classes of problems are considered. The first class pertains to 
flows wherein the basic (unperturbed) flow extends from x = -<» to 
x = -foo. The second class pertains to flows wherein there are two basic 
uniform flows separated by a discontinuity (such as a shock wave, contact 
surface, flame front, etc.) that moves with constant speed. The third 
class consists of several piecewise uniform basic flows, each separated 
by discontinuities moving with constant speed. For the latter class, 
attention is focused on cases wherein the discontinuities are centered. 
That is, they all originate at some fixed point, say x = 0, t = 0. Flow 
in an ideal shock tube (assuming the expansion wave has negligible thick- 
ness) is an example of a piecewise uniform flow with centered 
discontinuities . 


Basic Flow Infinite in Extent 


The solution of problems wherein the basic flow is infinite in ex- 
tent is given directly by equations (ll) and (12). Since the physical 
nature of this solution has already been discussed in detail, only two 
simple examples are treated herein. 


First, consider a flow for which AA = AA(0; m = q = f = 0. That 
is, the perturbations are due to a steady-state area variation in the 
tube. Integration of equations (lie) and (llf) gives (with 

|X = ^ AA.(-«») = AA.(+co) = 0) 


Ap + (x) r M AA(x) 

p = " 2 1 + M A 

Ap"( x ) _ r m aa(x) 

p 2 1 - M A 


(19) 


which define the acoustic waves at station x. The net pressure and ve- 
locity perturbations at this section are then 


Ap(x) = rM 2 AA(x) 

P 1 - M 2 A 

Au(x) -rM AA(x) 


( 20 ) 
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These are the same results as obtained directly from steady-state isen- 
tropic flow equations (e.g., ref. 5). Hence, steady flow through a 
slightly nonuniform tube can be viewed as the standing wave resulting 
from the superposition of unsteady acoustic waves generated by the source 

distribution p. = — df^‘ Note that the perturbations at station x 

depend only on the local area perturbation AA(x) (eqs. (20)). The ori- 
gin of these perturbations may be seen more clearly by considering a tube 
wherein the net area perturbation AA(x) occurs at station g A and an 
equal and opposite area perturbation occurs at £ B , as indicated in fig- 
ure 2(a). The equivalent planar source strength p at § A and | B is 

-puAA(x)/A and puAA(x)/A, respectively. For a subsonic basic flow, 
the pressure perturbations at (x,t) are due to the Ap + wave originating 
from the source at (£ a ,t a ) and the Ap" wave from the source at (| b ,t b ) 

(fig. 2(b)). For a supersonic basic flow, the pressure perturbation at 
(x,t) is due to the Ap + wave from the source at (i^T^) and the Ap" 

wave from the source at (| a ,t^) (fig. 2(c)). 


The second example is as follows i Coisider a gas to be stationary 
and uniform in a tube of constant area. Ab time r = 0, volumetric body 
forces of unit strength are distributed along the tube from g = 0 to 
£ = 1, and they remain at constant strength thereafter. That is. 


f(6,T) 

= 0 

for 

all 

i and 

f(e,T) 

= 0 

for 

i < 

o. g > i 


= 1 

for 

o < 

t < i 


(21a) 

(21b) 


The problem is to find the resulting perturbations. These momentum 
sources occupy the crosshatched region in figure 3(a). The perturba- 
tions at a typical point (x,t) arise from '.hose portions of the two 
characteristic lines through (x,t) which pass through the crosshatched 
region. (These portions are darkened in figure 3(a) for the two typical 
^°+^ S indicated therein. ) From equations (ll) it is seen that 

Ap (x, t) is positive and is proportional tc the length of the downstream 
propagating characteristic intersecting the crosshatched area. Also, 

Ap (x,t) is negative and is proportional tc the length of the upstream 
propagating characteristic intersecting the crosshatched area. As a re- 
sult, each of the numbered regions in figuie 3(b) has a different ex- 
pression for the local net pressure perturbation. The pressure distri— 
bution in the tube at times t = t * and t = t* 1 is indicated in fig- 
ures 3(c) and (d) . The solution is antisymmetric about x = l/2. This 
type of approach is applicable for finding perturbations induced by im- 
pulsive application of a magnetic field to a conducting field flowing 
through a uniform tube. 
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Basic Flow Containing a Single Discontinuity 


Consider a basic flow consisting of two regions, 1 and 2, separated 
by a discontinuity moving with velocity w (fig- 4(a)). Th e discontinu- 
ity may be a shock wave, contact surface, flame front, and so forth. The 
flow in regions 1 and 2 is assumed to be perturbed by mass, momentum, and 
heat sources. The problem is to find the perturbation at a typical point 
(x,t) . 


« 

The net perturbation at any point is found by summing the contribu- 
tion of all the elemental sources influencing the point. The sources 
contributing to a typical point in region 2 lie along the characteristics 
noted in figures 4(b) and (c) for the basic flows indicated therein. 

These characteristics are found in the following way. First, the down- 
stream and upstream propagating acoustic characteristics and the particle 
path characteristic are drawn through (x,t). These are lines a, b, and 
c, respectively. These characteristics intersect the discontinuity at 
points B and C. All the possible characteristic lines in regions 1 and 
2 which terminate at points B and C are then drawn. 0 This gives the ad- 
ditional characteristic lines used in figures 4(b) and (c). The physical 
significance of these lines is as follows. The sources along lines a, b, 
and c contribute directly to the perturbation at (x,t) in the manner dis- 
cussed in the previous sections. The sources along linetd generate a 
downstream propagating acoustic wave which arrives at (ig,^)- The lat- 
ter is designated Ap^g. This interacts with the discontinuity and 
generates a reflected wave (designated Ap 2 ^g) which propagates along 

line b so as to arrive at (x,t) . Similarly, the sources along lines 
f and j each generate a disturbance at point B which interacts with the 
discontinuity to contribute to the upstream propagating acoustic wave _ 
A P ; *. The sources along g, h, i, and k generate a disturbance at point 
C which interacts with the discontinuity so as to contribute to the en- 
tropy perturbation at point C in region 2 (As 2>c ), which then propagates 

along line C so as to arrive at (x,t). 


The ratio of the generated wave (at a discontinuity) to the incident 
wave is termed a "transfer function" herein. Acoustic reflection and 
transmission coefficients are special cases. The numerical value of the 
transfer function depends on the nature of the incident wave, the gener- 
ated wave, and the discontinuity. Transfer functions for various situa- 
tions are given in appendix C. For example, the following transfer func- 
tions must be known at point Bin figure 4(b) in order to compute the 
perturbation at (x,t) : 

• ‘ " ' ' The net perturbation Ap“ B /p ? is the sum 


and (Ap 2 yp 2 )/(AS 1 , B /Cv,i ) - 


3 P. 220 of ref. 12 gives all the possibilities. 
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of the contributions from the three incident waves noted in the denom- 
inators of the latter three transfer functions. 


Formulas for the characteristic lines and intersection points in 
figures 4(b) and (c) are given as follows. The coordinate system is 
chosen so that the discontinuity goes through the origin. 


Line a: 

! = 

x - (ug + ag) (t - t) 

x=t-[ 

(x - l)/(u 2 + ag)J 

(22a) 


== 

^ T ) 

s *a<&> 



Line b: 

i = 

x - (ug - ag)(t - t) 

T=t - [ 

(* - |)/(u 2 - a 2 )J 

(22b) 


s 

£ b ( T ) 

s ^bCs) 



Line c : 

t = 

X - Ug(t - t) 

T=t - [ 

(x - |)/Ugj 

(22c) 


s 

^c(t) 

s T c(0 



Line d: 

1 = 

i B - (u 2 + a 2 )(T B - t) 

T = T B - 

[Ub ■ £)/ ( u 2 + a 2) 

(22d) 


= 


= T d^) 



Line e : 

1 - 

- (u B - s 1 )(t b - t) 

T = T B - 

[(^B " 0/ (up - a 1 ) 

] (22e) 


= 

ft e (T) 

- T e (i) 



Line f: 

i = 

" u f T B " T ) 

x = r B - 

[d B - 0/upJ 

(22f) 


= 

if*( T ) 

- 



Line j : 

£ = 

E b - (up + a 1 )(x B - t) 

x = r B - 

[(£ B - 0/ (Up+ ap) 

|(22g) 


= 


h r jC |) 



Point B: 

^B = 

f x " ( u 2 - **2^1 

x - (ug - a 2 )t 

(22h) 


l_w - (u 2 - a 2 ) J w 

TB w - ( u 2 " a 2) ‘ 

Point C : 

» 

w(x - Ugt) 

X — Upt 

(221) 


w - Ug 

U V - 

u 2 

The equations 
lines d, e, f, 

for lines g, h, i, and 
and j by replacing the 

k are found 
subscript B 

from the equations for 
with the subscript C. 


The perturbations at (x,t), in figure ..-(b), can now be expressed 
formally. Use the following notation for the integrals appearing in 
equations (l2a) and (12b): 
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+ _ I5_ 

24 a " 2P 2 


r =Ii- 

^e" 2p l 


f ^ % <Lz(lS) 

f 2 (l»] 

|d 2 U^) + ( Cp T) 2 

a 2 Jfi-fi a (T) 

r , % ^s) 

fpU^n 

L p 1 (| ^ t) + (c p T) x ' 

“i Ji=t e (t) 


dx (23a) 


dr (23b) 


and so forth. The pressure perturbations at (x,t) can then be written 
(from eqs . (12) and fig. 4(b)) 


Ap^(x,t) 

P2 


Ap~(x,t) 


w 


/ \t Ap 2,li/ p 2 L* V B + 

Apt Jv, dL 


(24a) 


Ap 2,B^ P 2 

Api ?B / p l 


B ^2,B' ^2 

(^.leC + / t 4pU ' T) ] 5= 


dT 


The entropy perturbation at (x,t) is 

Asg(x,t) ^ As 2?c ^ 1 

c v,2 c v , 2 + ^ C v pT) 2 


J [a 2 (^ T )J| = | c 


(24b) 

dr (24c) 


■where 


As 2 ,c _ ^2 


'v,2 Ap+ J p 


A s 2 ,c/ c v j 2 A VC + A VC + 

2,(/ P 2 ' £ ’ 6e '-” Ap iV Pp ' ’ h ''-“ 


- 


As i,c/ C v ; i ^ C v pT d 


Sc 

J [aiU,T)] |= ^ 


dT 


The net pressure, velocity, density, and temperature perturbations at 
(x,t) can then be found from equations (11a) to (lid}. 
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. ® le ex P r essions for the perturbations at (x,t) in figure 4(c) re- 
quire two . additional terms because of charecteristic lines j and k. The 
perturbations at a point in region 1 are fcund in a similar fashion. 


As an example, consider the propagation of a shock wave in a duct 

having small area variations. 
This problem has been discussed, 
from other viewpoints, in ref- 
erences 6 to 8. The wave dia- 
gram for obtaining the pertur- 
bations behind the shock (region 
2) is indicated in sketch (g) . 

- fluid in region 1 is assumed 
to be uniform and at rest so 
that no disturbances originate 
in that region. Thus, with 
AA = AA(0, m = q = f = 0, and 
AA(-oo) = o, equations (24a) and 
(24b) give 



-TtM. 


2(1 + Mo) 


AA(x) 

A 


(25a) 


Ap 2 (x,t) 
P 2 


= I 


A P2,b/p2 

Ap 2,b/ P 2 


2U a ). 


B 


12 M 2 

w 




| ~AA(x) _ ^(6 B ) 1 A P2,B /p 2 r 2^2 AA(| b ) 

*- A A J Ap^g/pg L 2 d + M ?) A 


(25b) 


where (Ap^g/p^/tAp^g/pg) is found from equation (ClOa) . The entropy 
perturbation at (x,t) is found by evaluating equation (25a) at (L,r-,). 
which gives Ap 2 ^ c /p 2 , and multiplying by th ; transfer function 
(As^ c /Cv^ 2 )/(Ap^ c /p 2 ), as given by equatio i (ClOb). The present solu- 
tion is in agreement with the previous- treatments of this problem. 


Basic Flow Containing Centered Discontinuities 

Basic flows which contain a number of discontinuities, but are uni- 
form between these discontinuities, can be treated by further extending 
the methods of the previous sections. In the present section, examples" 
are considered wherein the discontinuities a:~e centered at x = 0, t = 0 
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The first example is that of a piston starting impulsively from 
rest and moving with uniform velocity -Up thereafter (fig. 5). The 
fluid in the tube is initially at rest (region l). A shock is generated 
such that the fluid between the piston and the shock moves with a basic 
velocity equal to that of the piston (i.e. , u 2 = Up). The fluid in 
region 2 is assumed to be perturbed by mass, momentum, and heat sources, 
and the resulting perturbations are desired. Wo sources are assumed in 
region 1. The characteristics influencing a typical point in region 2 
are indicated in figure 5(b) . There are an infinite number of charac- 
teristic line segments due to successive wave reflections at the shock 
and piston. The expressions for the pressure perturbation at (x,t) are 
of the form: 



(26b) 

The entropy perturbation can be found by evaluating equation (26a) at 
point C (so as to have Ap+ Q ) and then utilizing equation (ClOb). 

4 

Since the piston moves with constant velocity, the reflection co- 
efficient at the piston ((Apg^/pgV^p-j^A/^) ’ ^P2 ,d/p2^/(AP2,d/p2^ > 
etc.) equals 1 (eq. (C5)). The reflection coefficient at the shock 

((AP2 ?b / p 2^/^ Ap 2,b/ P 2^ ^ Ap 2,e/ P 2^/^ Ap 2,e/ P 2 ^ ’ etc ') ls zero for 

M = u /a-, = 1 (eq. (ClOa) ) . For T near 1.4, it is small for all 
values of M . In the latter cases, only the first few terms are re- 
quired in equations (26) * Figure 5(c) indicates the characteristic 
lines which are considered when the reflection coefficient at the shock 
is essentially zero. 

4jf the piston velocity has small nonuniformities AUp, the problem 
can be treated as the superposition of the constant Up case plus the 
case wherein pressure waves of magnitude a P2,a/ p 2 = pAu p , a/ a 2 ’ 

Ap 2 r/ p 2 = TAu p r/ a 2’ etC ' are generated at points A, D, F, . . .. The 
latter follow from the acoustic relations. 
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The piston-driven shock problem was treated in reference 6, from 
another viewpoint, to obtain the effect of small variations of piston 
velocity and the effect of a small linear area variation. The effect 
of an arbitrary area variation (e.g., wall boundary-layer effect), as 
well as arbitrary heat and momentum addition, can be treated by the 
present method. 


The second example is the determination of nonuniformities in shock 
tubes due to unsteady wall boundary-layer action. This problem was 
treaoed in references 2 and 3, and the solution is summarized herein. 

Let regions 1 and 4 be the low- and high-pressure sections of a shock 
tube (fig. 6(a)). When the diaphragm breaks, a shock wave moves into 
region 1 with velocity u g while an expansion wave (assumed to have 
negligible thickness) propagates into regicn 4 with velocity -a 4 . A 
contact surface separates the shock compressed gas (region 2) from the 
expanded gas (region 3). Because of the fluid motion relative to the~ 
wall, a boundary layer develops along the wall between the expansion 
wave and the shock as indicated in figure 6(b). The problem is to find 
the perturbations (due to the boundary layer) at a typical point in 
region ^ . The characteristics considered in references 2 and 3 are in- 
dicated in figure 6(c). The pressure perturbation at (x,t) is then found 
from (in the present notation) 


Ap^(x,t) 


Po 


; V + A P2,a/p 2 / _ y A + AP2 ,a/p 2 

a P2,a/? 2 V 2 ^e/T E a P3^ a /p 3 


• 3,1 


d/T 


D 


(27a) 


Ap 0 (x,t) 


Pn 


= I 


1 + A P2,B/P2 L+ \ t B + A P2,F/P2 

^b/T B Ap 2 ^ B /p 2 _\ 2 ^f/ TF Ap 2 ^ F /p 2 


2,i 


h/Tn 


H 


t p' 1 


Ap 2,f/p2 

AP 3,F /P 2 


(27b) 


Tne bracketed term on the right side of equation (27b) equals A Pg n / P a 
and may be found by evaluating equation (27a) at point B. The entropy 
perturbation at (x,t) is found by evaluating equation (27a) at point C 
and applying equation (ClOb). The source distribution in the integrals 
of equations (27) is given by 6 


P 2 - p 2 Zv 2 /A 

^3 = P3 Z V A 


( 28 ) 
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where v 2 and are the normal velocities (positive when directed 

inward) at the outer edge of the boundary layer in regions 2 and 3, re- 
spectively (as discussed in connection with eq. (2)). These velocities 
have the form 


(la) = 

k 2 

(29a) 

(U S T - l )” 2 

Ua) = 

k 3 

(29b) 

(a 4 t + |) 5 



where K 2 and K 3 are constants and N = l/2 or l/5 for wholly laminar 
or wholly turbulent boundary layers, respectively (refs. 2, 9, and 10). 5 
By using these forms , the integrals in equations (27a) and (27b) are 
readily evaluated. A discussion of the resulting shock-tube nonuniform- 
ities is given in references 2 and 3. 

The linearized characteristic method developed in reference 4 is an 
alternate method for determining shock-tube nonuniformities (see discus- 
sion following eqs. (18)). However, the source distribution in reference 
4 was based, in effect, on f(|,x) (obtained by averaging the boundary- 
layer wall shear across the tube cross section) and q(£,t) (obtained by 
averaging the heat transfer at the wall and the dissipation within the 
boundary layer) as opposed to the use herein of p(£,t) (i.e., effective 
area change due to boundary layer). Since f(£,T) introduces an anti- 
symmetric wave pattern (with respect to sign) , the results of reference 
4 differ by more than just a factor of proportionality from the results 
of references 2 and 3 (particularly with regard to the perturbations near 
the contact surface). This point is discussed further in references 2 
and 3. For the case of relatively thin wall boundary layers, it is clear 
that the source distribution used in references 2 and 3 is the correct 
one. The viewpoint of reference 4 may have some validity, however, when 
the boundary layer spans the entire tube cross section (provided the 
proper boundary-layer theory is used to obtain the wall shear, heat 
transfer, and. dissipation and provided the averaged perturbations are 
sufficiently small to justify a linearized approach) . The authors of 
reference 4 have extended their work, and the results are given in ref- 
erence 11. 

The flow resulting when a weak deflagration is initiated at the 
closed end of a tube is another example of a basic flow containing 


'The expression for v-,(£,t) is somewhat in error because of neglect 
of the f;! n.L to width of the expansion wave and because of the presence of 
gas, originally from region 2, at the wall between the diaphragm location 
( x -- C.i) and she con too i: surface (see ref. 10 ) . Use of eq. (29b) seems 
adequate dor vest ujanoses (except, possibly, for strong shocks). 
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centered discontinuities (see p. 225 of re..’. 12). The basic flow con- 
sists of a shock wave, followed by a planar flame front, both propagat- 
ing with constant speed (fig. 7). The shock induces a velocity u 2 
in region 2. The closed-end boundary condition requires u 3 = 0. (For 
a given flame, the shock strength is determined by the boundary condi- 
tion u 3 = 0.) The effect of the wall boundary layer on the flame and 
shock propagation can be treated in a manner similar to the shock-tube 
problem. Thus, the source distribution is taken to be p = plv/A, where 
v is found from boundary-layer theory. Note that the boundary layer in 
region 2 is the same as that in region 2 of the shock tube. However, 
since Uj =0 in the present example, the boundary layer in region 3 is 
essentially a thermal boundary layer only. The solution of this thermal 
boundary layer can be found using the methods of reference 9 provided 
the flow across the flame front is assumed uniform (i.e., the boundary 
layer from region 2 is neglected) . The characteristics contributing to 
the perturbations at (x,t) in region 2 are indicated in figure 7(c). It 
is assumed therein that the reflection coefficient at the shock is zero 
(appropriate for weak shocks). Transfer fi notions, for use at a planar 
flame front, are derived in appendix C. The effect of flame front dis- 
tortion and the corresponding changes in the mean flame speed are not 
considered in the present formulation. The effect of temperature and 
turbulence level on flame speed is also neglected. 


CONCLUDING REMARKS 

A source distribution method has been presented for determining 
perturbations due to small unsteady area variations, mass, momentum, and 
heat additions in a uniform (or piecewise uniform) basic one— dimensional 
flow . 


In the present method, the perturbed flow field is decomposed into 
three mutually independent wave systems, namely (1) upstream propagating 
acoustic waves Ap~, (2) downstream propagsting acoustic waves Ap + , and 
(o) entropy waves As . The pressure and velocity perturbations at a 
given section depend only on the local values of Ap“ and Ap + and are 
independent of the local value of As. However, a knowledge of As is 
required if the local density and temperature perturbations are to be 
found. The present method is compared with the linearized method of 
characteristics (e.g., ref. 4) in the paragraph following equations (18). 
The linearized method of characteristics employs AP, AQ, and As as 
Uie dependent variables. But AP and AQ are not independent of As 
for nonisentropic flows. Hence, AP, AQ, ani As do not form a system 
of three mutually independent waves as do Ap + , Ap~, and As. 



Several alternate viewpoints can be used to solve the problems con- 
sidered in the body of the report (e.g., refs. 4, 6, 7, and 8). It is 
hoped that the present elementary discussion adds additional physical 
insight for these flows. 


Lewis Research Center 

Rational Aeronautics and Space Administration 
Cleveland, Ohio, February 6, 1959 
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APPENDIX A 
SYMBOLS 

cross-sectional area of tube 
points on characteristic lines 
speed of sound 
characteristic lines 
specific heat at constant pressure 
specific heat at constant volume 

body force, per unit volume, acting in +x direction 
see eqs . (23a) and (23b) 
perimeter of tube cross section 
Mach number, u/a 

mass addition, per unit volune, per unit time 
Riemann variables (eqs. (13d) and (l3e)) 
pressure 

heat addition, per unit volune, per unit time 

gas constant 

entropy 

temperature 

time 

velocity in +x direction 

velocity (normal to wall) of fluid particle at wall, or 
at outer edge of thin boundary layer (positive inward) 

velocity of discontinuity 
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X 

distance along tube, taken in stream direction 

r 

ratio of specific heats 

A 

prefix denoting perturbation 

A 

( c p T W ( c p T b 

M- 

equivalent mass addition, per unit volume, per unit time 


distance along tube, defining source point 

P 

density 

T 

time, defining source point 

Subscripts : 


1,2,3,.. - 

regions 

A jB y C y . . • 

points on characteristic lines 

f 

flame front 

s 

shock wave 

Superscripts : 


( ) + 

perturbation associated with acoustic wave propagating in 
+x direction 

( r 

perturbation associated with acoustic wave propagating in 
-x direction 

D 

used with yx, q, and f to indicate source strength defined 
on a per unit cross-sectional area basis 

Special Notation: 

B±( ) 

Bt 

eq. (I3f) 

D( ) 
Dt 

eq. (I3g) 
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Example : 



eqs. (23a) and (23b) 


pressure perturbation at point A, 
with acoustic "wave propagating 


region 2, associated 
in +x direction 
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APPENDIX B 


ELEMENTAL SOURCE SOLUTION FOR M > 1 

The wave due to an elemental source in a subsonic stream was derived 
in the body of the report. The case of an elemental source in a super- 
sonic stream is treated herein. 


Planar mass , momentum, and heat sources are placed at x = 0 


t = 0 and remain at 

/ r. 



(h) 


at 

constant strength thereafter. The resulting t-x 

diagram, for M > 1, is shown in sketch 
(h) . Line a is a downstream propagat- 
ing acoustic wave, line b is a contact 
surface, line c is an upstream propagat- 
ing acoustic wave (which is swept down- 
stream, since M > l), and line x = 0 
is a line across which discontinuities 
also originate. The perturbations are 
constant in each of regions 2 to 4. 

These are found by considering the per- 
turbations arising from each of the 
above-mentioned lines . 


Line a : The acoustic relations give 


Line b : 


Line c ; 


Ap^ = paAu 2 I 

Apg = — Au 9 l 
a ^ I 

Since this line is a contact surface, 
Ap 3 = Ap 2 I 

Au^ = Aug j 

From acoustic relations. 


Ap 3 - Ap 4 = -pa(Au 3 - Au 4 ) 
Ap 3 “ a ~ ■ ~ ~ (Au 3 — *Au.,< ) 


(Bl) 


(B2) 


(B3) 
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Line x = 0 : Integration of equation;; (l) across x = 0 gives 

uAp 4 + pAu 4 = p 

uAu 4 + (l/p)Ap 4 = f/p 


u 


rA Pj _ ^P4\ = q 

iP ~ T P I ’' V P T 


(B4) 


Equations (Bl) to (B4) are* nine equations in nine unknowns. The solu- 
tions can be written 


4P2 _ ^3 _ r 

p ~~ p ~~ 2pa(l + M) 
Ap 4 


v 1 + c p t + V 


ym 


P pa(jVT - 1) 

_L_ 


p + 


CpT Ma 1 


Asj As 4 


c v c v 


c v puT 


(B5) 


The other perturbations can then be 



found from equations (Bl) to (B4) . 

New assume that equal and oppo- 
site (with respect to the previous 
example ) planar mass momentum and 
heat sc urces are placed at x = 0 
at t = At and remain at constant 
strenglh thereafter. The t-x dia- 
gram is given in sketch (i). The 
lines *re displaced upward by an 
amount At from those in sketch (h) 
The perturbations in regions 2 { , 3 T , 
and 4 r are equal and opposite to 
those jn regions 2, 3, and 4. 


( ' ) 

The flow field associated with an elemental mass, momentum, and 
heat source at x = 0, t = 0 is found by 3 inearly superposing the flows 
in sketches (h) and (i). The superposition is indicated in sketch (j) . 
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r 4 + 4‘ 



The perturbations are zero everywhere except between lines a 
and c-c’ (since the perturbations in regions 2, 3, and 4 are 
opposite to those in regions 2', 3 ' , and 4', respectively), 
bations between lines a-a' can be written 


Between lines c-c' 7 


Ap H 


T 


2pa(l + M) 


c T + 


Au‘ t 


Apr 

p 


1 Ap] 

r p 


the net perturbations are 


Ap" 


r 


2pa(M - 1) 


_ , _JL_ 

^ c T 


Au~ 

a 


Ap' 


1 

r p 




Between lines b-b', the net perturbations are Ap = Au = 0 

As _ q 
c v = c v puT 


Ap 

P 


As 

c 


-a’ , b-b T , 
equal and 
The pertur- 


(B6a) 


(B6b) 


and 


(B6c) 
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By considering the limit as At -*■ 0, equations (B6) define the perturba- 
tions due to an elemental planar source at x = 0, t = 0 in a supersonic 
stream. These equations differ from those for the M < 1 case (i.e., 
eqs. (9)) only in the terms (l - M), (M - l) appearing in the expressions 
for Ap . Use of jl - M| in equations (10 makes the latter applicable 
for both M < 1 and M > 1 . 

It may be of interest to interpret the flow in sketch (h) in terms 
of the elemental source solutions. For example, the net pressure per- 
turbation at an arbitrary point in region 4 (see sketch (k) ) may be in- 
terpreted as consisting of 
Zip = Ap + + Ap“ where Ap + is 
the downstream acoustic wave 
c riginating from the source at 
! = 0, t'" = t -[x/(u+ a)] 
while Ap - is the upstream 
propagating acoustic wave orig- 
inating from the source at 
l = 0, t' = t -[x/(u - a)]. 

The entropy perturbation at 
(x,t) is equal to the entropy 
perturbation which originated 
irom the source at | = 0, 

X, £ i" = t - (x/u). The greater 
generality of the elemental 
source viewpoint is illustrated 
by considering the source strength at £ = 0 to vary with time. Then, 
the perturbations at any point can be founc directly from equations (B6) 
provided the right sides are evaluated at t'", t 1 , and x", respectively. 
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APPENDIX C 


TRANSFER FUNCTIONS 

In the body of the report It was necessary to know the ratio of the 
generated wave to the incident wave when an acoustic wave, or entropy 
wave, impinged on a discontinuity. Such ratios were termed "transfer 
functions." The transfer functions for various discontinuities and in- 
cident waves are found herein. 


Contact Surface 

The requirement that pressure and velocity be continuous across a 
contact surface leads to (e.g., ref. 2): 


Incident wave Ap^ : 



Incident wave Apg i 
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Constant-Velocity Discontinuity 

In this case the velocity of the discontinuity is assumed to be 
unperturbed. A piston moving at constant velocity is an example. 



Constant-Pressure Discontinuity 


In this case, the pressure at the discontinuity is assumed to be 
unperturbed. The open end of a tube is an example. 



A 4/ p 2 


(C6) 


The derivation of the transfer functicns for this case will be out- 
lined. Let M s = (u - u l)/ a l designate the shock Mach number relative 
to the fluid in region 1. 



u 2 



u , 


(p) 
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Normal shock relations give 

P2 _ ^s - (r - i) 
Pi " r + i 


Uo - Ui 


2 M s - 1 


r + i m c 


s 2 - S 1 


In 


’ 

2rMg - (r - i) 

1 

1 

co m 

i—l 

+ 

I 

r 


x + 1 

1 

j_(r - i)Mg + 2 



(C7a) 


(C7h) 


(C7c) 


Taking the differential of equations (C7a) and (C7h), eliminating 

Aa i 1 f r - 1 Ap i 1 


and noting 


an 


2 y r pi r 


?) 


yield 


Ap 2 

Ap x 


f2 

/ Au 2\ 

Aui r - 1 M s ' 1 / Ap l 1 

P2 

" Pi 

z 

a l 

l r ^r) 

- r - r + i m s \ Pl + r - i c v j 


(C8) 


where 


2(r + i)m: 


Z = 


[Irid - (r - i)] + i) 

The differential of equation (C7c) , upon elimination of AM g , yields 

4 


As g - As^ ^ APg Ap-^ 

" pi 


) 


(C9) 


where 


(r 

- l) (M 2 - l) 2 

M 2 

s 

(r - Dm 2 + 2 ] 


Use of the acoustic relations rAu + /a = Ap + /p and IfAu - / a = -Ap“/p in 
equations (C8) and (C9) then yields the following transfer functions: 
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(r) 




37 


Incident wave As,/c ir : 



M? - 1 


Ap"/p- 

As^/c v = " (r + i')m s 


1 + z 


a l 


(C13a) 


As 2 /c v 

^l/^v 


= 1 + 


WM, 

ymD v j ; 


(C13b) 


PI anar -Flame -Front Discontinuity 


Transfer functions for a planar flame front are derived herein by 
assuming : 



(2) The flame speed relative to the gas in region 1 , u f " u ]_j -** s 
constant 

(3) The heat release at the flame per unit mass, per unit cross- 
sectional area (denoted Q ) , is constant. 


The e on it ions 0 f 

motion, relative 

to the flame . are 


Continuity : 

Pi Of “ ,,J l) = Pp 

( u [* - u 2 ) 

(ci:m) 

Mormn turn : 

pi =■• Pu 


(Cl 4b) 

Energy : 

(% T >i + 5 = 

(y T) 2 

(C14c) 

g L • ci. O e m 

Pl = P^iTj I 

= p t: T : 

(C14d) 
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Assumption (l) is incorporated in equation; (C14b) and (C14c). The dif- 


' u f - 

u.j) ~ 0 or 

Au^ = 

: Au^, ind 

(3 ) > i. 

App 

Ap 2 ap 

__ _1_ 

(r Au - 

a 2 r l 


p 2 u f " 

u 2 n 

V 1 ®i 

¥z r! 



App 

Ap 2 



ATp AT 2 


^2 ^1 /A - l\ ^1 


fX - l\ ^1 

V A ) T i 


(Clba) 

(C15b) 

(C15c ) 
(C15d) 


where A = ( c pT) 2 /( c pT)p. Eliminating the density terms between equa- 
tions (C15a) and (C15d) yields 


r l Z ^1 a 2 r l Au 2 _ 

Yp — 1 Tp ap x 2 a., ^"l a^ 


(C16) 


where 


Z = (r 2 - x) (A - Dp-^) 

But ATp/Tp = ^(rp - l)/r 1 l(App/p 1 ) + (l/r] ) (Asp/c v ^p) ^ so that equation 
(C16) can he written 


i ^i \ a 2n 


Pl + rq - 1 c vA J apr 2 r 2 a 2 “ r l ap 


Also, from the definition of entropy, 

Asg _ Ap 2 Ap 2 ^ ^ Ap 2 

c v,2 p 2 TZ p 2 T p 2 


(C17) 


2 ' P 2 + T2 


= (1 - To) 


Api r 2 ATp 


App p T 2 ^1 

= Y + T 

Pl A rp c Vfl 


(C18) 
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where 


Y = 


(r x - i) 


1 ^2 


- 1)" 

ATi " 


- 1)J 


Equations (C15h) , (C17), and (C18) together with the acoustic relations 
then yield the following transfer functions : 


Incident wave Apt : 

s* i 



Incident wave Ap ^ • 



+ / 1 
Ap£/ Pq 

Ti a 2 

(C19a) 

Ap i /p i i 

Ti a 2 

+ z> + „ _ 

^2 a l 

A P2/l’2 

Ap^/Pq 

A Pl/Pl 
+ Ap£/p x 

(C19b) 

As 2/ c v,2 

Ap£/P 2 

Y 

(C19c) 

Ap£/p-l 

Ap£/ Pq 

Ap£/p 2 

n a 2 

1 + Z - rr- ~ 

r 2 a i 

(G20a) 

Ap+/p 2 

r i a 2 

1 + Z + ^ — 

r 2 a i 

Apt/ P-, 

— = 1 
Ap 2 /p 2 

Ap£/ p 2 

Ap 2 /p 2 

(C20b) 

As 2/ c v,2 

A Pl/Pl 

Y 

(C20c ) 

Ap 2 /p 2 

Ap 2 / p 2 
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(a) Area variation. 

(b) Sources affecting point (x,t) for subsonic flow. 

(c) Sources affecting point (x,t) for supersonic flow. 


Figure 2. - Source distribution method for finding steady 




(a) Basic flow containing discontinuity. 

(b) Characteristics contributing to perturbations at (x,t) for 
case 0 < Mp < 1, 0 < M 2 < lj u^_ < w < u-^ + Ug < w. Since 
(w - u- L )/a 1 < 1, discontinuity moves with subsonic velocity 
relative to gas in region 1. 

Figure 4. - Characteristics for finding perturbations in flows 
containing a discontinuity. 





(£ C ,r C ) 


Discontinuity 




(c) Characteristics contributing to perturbations at (x,t) for 
case 0 < M]_ < 1, 0 < M 2 < 1, u x + a x < w, u g < w. Since 

(w - u 1 )/a 1 > 1 , discontinuity moves with supersonic velocity 
relative to gas in region 1. 

Figure 4. - Concluded. Characteristics for finding perturbations 
in flows containing a discontinuity. 






48 


Diaphragm 


(o) 




rWall boundary layer 


TZZZZZZZZZ^. ’MZZZZZ 


^^5Z2222 hJXm 


(b) 




(a) Shock tube before diaphragm rupture. 

(b ) Flow at time t = t 1 . 

(C L Gi ;: r + a f eriStiCS maMn ^ ma j° r contribution to perturbations 


Figure 6. - Calculation of nonuni formitie; and attenuation 
thickne- s tUbeS ' ExpanSl ° n Wave assurae <l to have negligible 
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Flome fronts /-Shock wave 



(b) 



(a) Basic flow (t = t 1 )- 

(b) Wall boundary layers. 

(c) Characteristics for determining perturbations at (x,t). 
Reflection coefficient at shock assumed negligible. 

Figure 7. - Determination of perturbations due to wall boundary 
layer when planar deflagration is initiated at closed end. 
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